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ABSTRACT 


The  most  general  case  of  plane  wave  propaga¬ 
tion  when  normal  and  shear  stresses  occur  simul¬ 
taneously,  is  considered  in  a  material  obeying  the 
V.  Mises  yield  condition. 

The  resulting  non-linear  differential  equa¬ 
tions  have  not  previously  been  solved  for  any 
boundary  value  problem,  except  for  special  situa¬ 
tions  where  the  differential  equations  degenerate 
into  linear  ones.  In  the  present  paper,  the 
stresses  in  a  half-space,  due  to  a  uniformly  dis¬ 
tributed  step  load  of  pressure  and  shear  on  the 
surface,  are  obtained  in  closed  form. 
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INTRODUCTION 


I . 


The  basic  differential  equations  governing  wave  propa¬ 
gation  in  an  ideal  elastic-plastic  medium  subject  to  the 
V.  Mises  yield  condition  are  easily  formulated.  It  is  es¬ 
tablished  [Ij*  that  these  equations  are  hyperbolic  in  char¬ 
acter  ^  and  that  there  are  two  non-vanishing  pairs  of  charac¬ 
teristic  velocities,  +V  ,  +V  ,  However,  the  actual  solu- 

-  1  -  2 

tion  of  boundary  value  problems  is  complicated  by  two  causes 
in  general,  the  differential  equations  are  non-linear,  so 
that  the  are  functions  of  the  stresses,  and  there  are, 

a  priori  unknown,  moving  boundaries  between  elastic  and 
plastic  regions.  As  may  be  seen  from  a  review  of  the  field, 
[2J,  the  published  literature  contains  solutions  only  to  one 
dimensional  problems  of  plane,  spherical  or  cylindrical  pres 
sure  waves  where  just  one  of  the  two  characteristic  veloci¬ 
ties  enters  the  solution,  and  where,  furthermore,  the  dif¬ 
ferential  equations  degenerate  into  linear  ones,  so  that 
this  velocity,  +V,  is  a  constant.  No  multi-dimensional  prob 
lem,  necessarily  involving  non-linear  differential  equations 
with  two  pairs  of  non-constant  characteristics,  has  at  pres¬ 
ent  been  treated. 

As  a  step  towards  multi-dimensional  problems  the  pres¬ 
ent  paper  considers  a  case  in  which,  while  one-dimensional, 
the  differential  equations  are  nonlinear  and  both  character¬ 
istic  velocities  enter  the  general  solution.  Subjecting  a 
half-space.  Fig.  1,  to  a  suddenly  applied,  uniformly  dis¬ 
tributed  normal  stress  (^^(t)  combined  with  a  shear  stress 

*  Numbers  in  parentheses  [ ]  refer  to  the  bibliography. 
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,  plane  waves  of  pressure  combined  with  shear  are  gen 
erated.  For  the  case  of  step  loads,  a  (t)  =ah(t),T  (t)  - 
^  closed  form  solutions  for  the  stresses  and  veloci¬ 
ties  can  be  obtained.  For  such  loads  dimensional  considera 
tions  make  it  possible  to  convert  the  partial  differential 
equations  into  a  set  of  simultaneous  ordinary  differential 
equations  which  can  be  solved  by  quadratures  in  terms  of 
elliptic  integrals  and/or  simpler  transcendental  functions. 
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II.  FORMULATION  OF  THE  BASIC  EQUATIONS 


The  basic  equations  governing  wave  propagation  in  an 
ideal  elastic-plastic  continuum  consist  of  a  yield  condition 
defining  the  plastic  state,  a  set  of  const ituitive  equations 
relating  stress,  strain  and/or  their  derivatives,  and  the 
equations  of  motion. 

Starting  with  the  yield  condition,  this  report  will 
consider  a  material  obeying  the  v.  Mises  condition 


f  = 


J' 

2 


-  =  0 


(2-1) 


where 


J‘ 


(2-2) 


is  the  second  invariant  of  the  stress  deviator,  and  k  /  0 
is  the  yield  stress  in  simple  shear. 


eij 


For  an  ideal  elastic-plastic  material  the  total  strain 

g 

is  composed  of  an  elastic  part  e . .  and  a  plastic  part 


ef.  ,  so  that  the  total  strain,  and  strain  rate  are,  re- 

ij  ^ 

spectively 


+  €' 


ij 


(2-3) 


(2-4) 


Employing  conventional  symbols,  listed  on  pages  iv  and 
V,  the  elastic  portion  of  the  strain  rate  is 
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(2-5) 


-1_ 

2G 


s  .  .  + 
T-J 


1  -  2v 
6g(i+v) 


j 


1 


while  the  plastic  portion,  obtained  from  the  concept  of  a 
plastic  potential  [5]  is 


(2-6) 


The  quantity  A  is  a  function  of  space  and  time  defining 
the  rate  of  energy  dissipation  at  the  respective  point  in 
space  and  time  and  is  necessarily  non-negative,  A  >  0  .  In 
regions  of  plastic  deformation  inherently  accompanied  by 
energy  dissipation,  A  must  be  positive 

A  >  0  .  (2-7) 

For  the  v.  Mises  yield  condition,  Eq.  (2-1), 


so  that 


s 


ij 


1 

2 


(u.  .+u. 

i,j  j 


)  =  ^ 
i  ’  2C 


+ 


1  -  2v 


2G  ij  6g(1+v) 


J  6  .  . 
1  ij 


+  As 


ij 


(2-8) 


where  the  strain  rates  are  expressed  in  terms  of  the  com¬ 
ponents  u^  of  the  velocities 


1  /• 
=  2 


1.  J 


+  u.  .) 


J.l' 


(2-9) 


Finally,  retaining  only  linear  terms  in  the  equations  of 
motion,  the  latter  become 


a 


Si 

p  It 


(2-10) 
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The  equations  stated  above ^  i.e.,  Eqs .  (2-1),  (2-8), 
and  (2-10),  apply  in  regions  of  plastic  deformation  accom¬ 
panied  by  energy  dissipation.  In  general,  there  will  also 
be  regions  where  at  the  particular  time  considered  no  plasti.c 
deformation  or  energy  dissipation  occurs.  In  such  regions, 
while  the  strain  rate-velocity  relations  (2-9)  and  the  equa¬ 
tions  of  motion  (2-10)  still  hold,  the  yield  condition  (2-1) 
and  the  stress-strain  relation  (2-8)  are  no  longer  valid. 
Instead,  in  "non-dissipative"  regions* the  total  strain  rate 
is 


while 


(2-11) 


Jj  <  k®  .  (2-12) 

The  general  equations  are  now  applied  to  the  following 
plane  strain  problem  in  Cartesian  coordinates  x,  y,  z.  The 
originally  stress-free  half-space  x  )>  0  ,  Fig.  1,  is  sub¬ 
jected  to  time  dependent,  but  uniform,  loads  on  its  surface 
X  =  0  , 


(0,y,z,t)  =  o^{t)" 

'rj^y(0,y,z,t)  =  T^(t)  > 

=  0  ^ 


(2-13) 


where  the  prescribed  loads  and  vanish  for  t  •(  0 . 
Because  of  the  nature  of  the  yield  condition  (2-1),  the  in¬ 
tensity  of  the  shear  is  subjected  to  the  restriction 


I  I  <  •  (2-13a) 

*  "Non-dissipative"  regions  in  the  sense  used  here  include 
neutral  plastic  regions  where  . 
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In  a  plane  strain  problem  the  stresses  t  and  t 
and  the  displacements  in  the  z-direction  vanish,  while  all 
other  quantities  are  independent  of  the  coordinate  z  .  In 
the  present  problem,  the  prescribed  loads  are  also  indepen¬ 
dent  of  y  ,  so  that  the  response  too  will  be  independent  of 
y  ,  and  all  quantities  will  be  functions  solely  of  x  and 
t.  When  writing  the  various  equations  in  scalar  form,  all 
derivatives  with  respect  to  y  or  z  will  therefore  vanish. 

In  regions  of  plastic  deformation  and  energy  dissipa¬ 
tion,  A  )>  0  ,  Eq.  (2-8)  gives  four  meaningful  scalar  equa¬ 
tions  : 


1  1  ,  1  -  2v  1  .  T.  o 

2G  X  6g(1+v)  *^1  X  dx 


55  \  1-  A  =  0 


k  "z  +  A  =  0 


2G 


T  +  At 


1.  ^ 
2  dx 


(2-l4a) 

(2-l4b) 

(2-l4c) 

(2-l4d) 


where  t  =  t^^  and  u  ,  v  are  the  x-  and  y-components 
of  the  velocity. 

The  equations  of  motion,  Fq.  (2-10),  after  the  substi¬ 
tution  a  =  s  +  :r  J  ,  give  the  two  non-trivial  relations 
X  X  3  1  ■ 


ds 

X 

dx 


dj 


du 


dx 


~  dt 


dr  dv 

S?  “  P  5t 


(2-15) 
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Subtracting  Eq.  (2-l4c)  from  Eq.  (2-l4b),  one  finds 


X 

2G 


+  A(s  -s  )  =  0  . 
y  "z. 


Initially,  at  t  =  0  ,  and  are  both  equal  to  zero, 

so  that  the  above  equation  requires  s  =  s  for  all  values 

y  z  ^ 

of  X  and  t.  Noting  that  by  definition  s  +  s  +  s^  =  0, 

X  y  ^ 

one  finds 


s 

y 


(2-16) 


which  is  used  to  eliminate 


and 


y  z 

Thus,  the  number  of  stress  variables  is  reduced  to 

three,  namely,  s  ,  J  ,  and  t  .  Simultaneously,  the 

1 

number  of  independent  constituitive  relations  is  also  re¬ 
duced  to  three.  A  convenient  set,  obtained  by  obvious  manip¬ 
ulation  of  Eqs.  (2-l4)  is 


1 

• 

6u 

6G 

J 

1 

dx 

1 

s 

+  As 

_  2 

ou 

>  (2-17) 

2G 

X 

X 

3 

Ox 

1 

<7* 

+  At 

1 

dv 

2G 

1 

~  2 

ox  J 

where 


«  3 

K 


A  f  M  \ 
^  \  J 


2(l+v) 
1  -  2v 


/  /-> 
V 


Using  Eq,  (2-16)  leads  finally  to  the  yield  condition 
in  the  simple  form 

I  s2  +  .  (2-19) 
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This  relation,  together  with  the  five  equations  (L-I5)  and 
(2-17),  are  available  to  determine  the  six  unknowns  s  ,  J  , 

X  1 

T,  u  ,  V  and  A  >  0  in  regions  of  plastic  energ^  dissipa¬ 
tion  which  will  be  referred  to  as  "dissipative"  regions. 

There  will  also  be  complementary  "non-dissipative"  regions 
without  plastic  deformation  or  energy  dissipation.  In  such 
regions  the  equality  (2-1)  is  to  be  replaced  by  the  inequality 

(2-12).  As  previously  noted,  the  differential  equations  for 

P 

the  non-dissipative  regions  may  be  obtained  by  using 
instead  of  Eq.  (2-8)  in  conjunction  with  Eqs .  (2-9).  This 
is  equivalent  to  substituting  X  =  0  into  Eqs.  (2-l4)  giving 


A  a.  1  -  2v  ;  _  ^u"^ 

2G  X  6g(1+v)  “Sx 


X 

2G 


1  -  2v 
6g(i+v) 


j 

1 


X 

2G 


1  -  2v 
6g( i+v) 


j 

1 


0 


(2-20) 


1-  r  -  i  dv 

2G  “  2 


Noting  that  Eq.  (2-l6)  is  still  valid  in  non-dissipative 
regions,  and  recalling  the  definition  of  p  ,  Eqs.  (2-20) 
may  be  reduced  to 


1 

T 

du 

PG 

J 

1 

~  dx 

1 

2  du 

2G 

s 

X 

~  3  dx 

1 

• 

1  dv 

2G 

T 

~  2 
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Moreover,  the  inequality  (2-12)  becomes 


^  sj  +  .  (2-22) 

Subject  to  this  inequality,  the  three  equations  (2-21)  and 

the  tv;o  equations  of  motion  (2-15)  form  a  system  of  five 

simultaneous  differential  equations  for  the  five  functions 

s  ,  J  ,  T  ,  u  and  v  .  This  system,  of  course,  is  equi- 
^  1 

valent  to  the  conventional  equations  governing  elastic  wave 
propagation . 

The  complete  boundary  value  problem  to  be  solved  involves, 
therefore,  an  a  priori  unknown  number  of  dissipative  and  non- 
dissipative  regions  with  unknown  and  moving  boundaries.  The 
solution  in  each  region  is  to  be  obtained  from  the  appropriate 
set  of  equations  given  above,  and  from  matching  and  external 
boundary  conditions. 

The  solution  of  the  boundary  value  problem  formulated 

above  for  a  general  input  C!^(t)  ,  T^(t)  is  a  formidable 

problem,  even  if  approached  by  numerical  methods.  As  stated 

in  the  introduction,  the  basic  set  of  non-linear  differential 

equations  is  hyperbolic,  and  has  two  pairs  of  characteristic 

velocities,  +V  .  and  +V  .  The  use  of  the  standard  numeri- 
-  1  -  2 

cal  approach  based  upon  the  method  of  characteristics  is, 
however,  complicated  by  unknown  moving  boundaries  between 
the  regions  and  by  the  occurrence  of  discontinuities. 

An  alternative  numerical  approach  is  available  which 
eliminates  the  difficulties  caused  by  moving  boundaries  and 
by  discontinuities.  It  replaces  the  differential  equations 
by  finite  difference  ones  and  uses  artificial  viscosity  terms, 
as  suggested  in  [^],  in  order  to  obtain  some  measure  of  ap¬ 
proximation  to  the  discontinuities  which  occur  in  the  exact 
solution . 
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It  should  also  be  noted  that  for  either  of  the  two  nu¬ 
merical  approaches,  a  starting  solution  valid  for  a  short 
time  is  required.  The  usual  power  series  approach  used  to 
obtain  a  starting  solution  is  unsuccessful  if  o^(t)  or 
T^(t)  is  discontinuous  at  t  =  0  ,  because  the  point  x  =  0, 
t  =  0  then  becomes  a  singularity. 

The  remainder  of  this  paper*  will  not  concern  itself  with 
the  application  of  either  of  the  mentioned  numerical  approachf^s 
to  the  general  problem,  but  will  consider  only  the  much  simp¬ 
ler  situation  when  the  applied  surface  loads  are  step  func¬ 
tions  in  time.  For  this  case  solutions  in  closed  form  can 
be  obtained.  In  addition  to  giving  some  idea  of  the  nature 
of  the  response  in  the  general  case  of  time  dependent  load¬ 
ing,  the  results  give  a  starting  solution  for  a  numerical 
treatment  of  the  general  problem. 


*  Except  for  a  discussion  of  open  questions  concerning  the 
application  of  the  method  of  characteristics  in  Appendix  B. 


-10- 


III.  ^lALF  SPACE  SUBJECTED  TO  COMBINED  NORMAL 


AND  TANGENTIAL  STEP  LOADS 


Consider  the  case  where  the  surface  loads  are  step  func¬ 
tions  of  time 


o 

o 


(t) 


=  o^h(t) 


T 

O 


It) 


Toh(t) 


(3-1) 


Using  dimensional  considerations,  the  number  of  dimensional 
quantities  appearing  in  the  equations  of  Sec.  II  may  be  re¬ 
duced  to  four,  i.e.,  x,  t,  G  and  p  ,  by  replacing  the  un¬ 
known  variables  s  ,  J  ,  t,  u  and  v  by  non-dimensional 

X  1 


(3-2) 


respectively,  while  A  is  expressed  by  a  new  function  A  , 
defined  by 

A  =  ^  A  .  (3-3) 

The  four  dimensional  quantities  may  be  arranged  into  a  single 
non-dimensional  independent  variable 


U  = 


(3-4) 


All  other  possible  non-dimensional  combinations  of  x,  t,  G 
and  p  are  functions  of  U  ,  and  are  therefore  equivalent 
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to  Eq.  (3”^).  ThuSj  the  non-dimensional  quantities  (3~2)  and 
(3“3)  solely  functions  of  the  independent  non-dimensional 

variable^  u  j  and  of  the  non-dimensional  parameters  b  j 
OqA  and  . 

A,  DISSIPATIVE  REGIONS 

To  convert  the  partial  differential  equations  obtained 
in  the  previous  section  into  ordinary  differential  equations 
with  respect  to  U  ,  note  from  Eq.  (3~^) 


±  A? 

dx  t  VG  dU 


a 

at 


u 

t  dU 


(3-5) 


(3-6) 


J 


Using  Eq.  (3-3)  ^  condition  X  )>  0  which  must  be 
satisfied  in  dissipative  regions  becomes 


UA'  <  0 

where  '  denotes  differentiation  with  respect  to  U  . 
The  relations  (2-17)  assume  the  form 


13  k 


0 


u  fx  +  !x  UA'  +  ^ 
2  k  +  k  +3 


U  T' 


ii'  =  0  ) 


2  k  +  ^  UA'  +  2  v'  =  0 


(3-7) 


(3-8) 


and  the  equations  of  motion,  Eqs.  (2-15),  become 
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s'  T  J 
_2i  +  1  _i 
k  3  k 


U  /l/p  u'  =  0 


(3-9) 


^  +  u  v'  =  0 


J 


The  yield  condition  (2-19)  in  terms  of  s  /k  and  'r/k 


IS 


i(!x 

4Vk 


+ 


_T 

,k 


=  1 


(3-10) 


and  its  derivative  with  respect  to  U  is 


(3-11) 


Equations  ('3“8);  (3-9)  ^nd  (3-II)  form  the  following  set  of 
simultaneous,  linear,  homogeneous  equations  for  the  deriva¬ 
tives  of  the  unknown  functions: 


"  0 

u 

6 

1 

0 

0 

0 

—  s 

X 

k 

u 

s 

j' 

0 
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0 

0 

X 

—L 

2 

3 

k 
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1 

1 

3 

u 

0 

0 

0 

¥ 

0 

0 

u 

1 

T 

0 

2 

2 

k 

k 

0 

0 

0 

1 

U 

0 

3. 

4  k 

0 

0 

T 

k 

0 

0 

ua’ 

=  0 


(3-12) 
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If  the  determinant  of  the  system  (3-12)  does  not  vanish, 
s'=j'  =u'  =t'  =v'  =UA'  =0  ,  and  Eq.  (3"7)  being  vio- 
lated,  the  above  equations  inherently  do  not  apply.  On  the 
other  hand,  solutions  where  UA'  <(  0  may  exist  if  the  de¬ 
terminant  of  Eqs.  (3-12)  vanishes,  which  gives  the  condition 

^  sj  (i-u2)0-3u^)  -  t2u2(2|+p-3u2)  =  0  .  (3-13) 

This  condition,  while  necessary,  is  not  sufficient  to  ensure 
UA*  <(  0  ,  which  remains  to  be  proved  later.  Equation  (3-13) 
expresses  the  variable  U  in  terms  of  the  stresses  s  and 
T  .  However,  since  s  and  t  are  related  by  the  yield 
condition  (3-10),  either  may  be  eliminated  with  the  result 


and 


s 

^  4-  Ai 

f4+6-3u2) 

ic  - 1  V 

31 

ip+U^) 

(3-14) 


T 

k 


=  + 


p  + 


(3-15) 


The  fact  that  U,  s  ,  and  t  must  be  real  results  in 
restrictions  on  the  extent  of  dissipative  regions.  The  bounds 
on  potential  locations  of  dissipative  regions  depend  on 
whether  the  value  of  P  is  greater  than,  less  than,  or  equal 
to  three : * 


*  For  the  usual  range  of  Poisson*  s  ratio,  0  ^  v  ^  ,  the 
value  p  ,  Eq.  (2-l8),  varies  from  2  to  oo  . 
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For  >  3:  0  <  <  1  or  (3-l6) 

For  fi  <3:  0  <  U^  <  U^  <  or  1  <U^<U^<  \j^  (3-17) 

For  =  3:  0  <  <  (\JY  (3-l8) 

where  and  u  designate  the  end  points  of  a  dissipative 

region.  Plot.?  of  Eqs .  (3-14)  and  (3-15)  typical  values 

of  p  sre  shown  in  Figs.  2  and  3j  respectively.  As  required 
by  the  two  inequalities  in  each  of  Eqs.  (3-l6,17)j  there  are 
two  separate  branches  except  for  the  special  case  P  =  3 
when  the  branches  for  s  merge,  while  the  ones  for  t  join 
in  a  cusp  at  U  =  1. 

It  is  interesting  to  note  that  the  bounds  U  =  V ( 4+P)/3 
and  U  =  1  ,  respectively,  correspond  to  the  velocities  of 
P-waves  and  S~waves  in  an  elastic  medium.  The  bound  u=  V p/3 
will  be  seen  in  Subsection  B  to  correspond  to  the  velocity 
of  propagation  of  plastic  shock  fronts. 

In  a  dissipative  region  where  the  determinant  of  Eqs. 
(3“12)  vanishes,  the  variables  j'  ,  \/^  u'  ,  x'  ,  s/pG  v' 

and  Ua'  can  be  written  in  terms  of  s' 

X 


j'  =  ^  s' 

1  3U'^  -  P  X 

(3-19a) 

(3-19b) 

T'  =  -  s' 

H  \  xy  X 

(3-19c) 

^  (f)  < 

(3-l9d) 

UA'  u(4+e-3u^)  3. 

2(3u-'-p)s^  ®x  • 

(3-19e) 
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This  set  of  equations  becomes  inapplicable  for  U  =  V 3/3  , 
where  T  =  0  .  The  physical  reason  for  the  breakdown  will 
become  apparent  in  Subsection  B. 

It  is  now  necessary  to  investigate  whether  Eq,  v3“7).' 
a  prerequisite  for  dissipative  regions,  is  satisfied  every¬ 
where  in  the  ranges  given  by  Eqs.  (3“l6  to  l8) ,  Using  Eq. 
(3-14)  to  express  u  and  s'  in  terms  of  U  ,  Eq.  (3“19e) 
becomes 


3[u^  +  p  +  -fr-'/F(P+TT][u2  +  b - ] 


UA  = 


sO+u^)(p-3u®) 


<  0 


(3-20) 


or 


U2  +  j3 - —  Vp(|3+l) 


3  -  3U^ 


<  0 


(3"21) 


Eq.  (3“21)  is  valid  if  either 


>  -3r-7p(p+l)  -  p  and  >  |  (3-22) 

•Jl  ^ 

or 


^  yp(e411  -  p  and  U“  <  I  .  (3-23) 

■Jl  ^ 

In  the  pertinent  range,  p  )>  2  ,  the  following  inequalities 
apply: 


For  p  >  3:  1  <  --^Vp(3+1)  -  P  <  I  (3-24) 

n/3  ^ 

For  p  <  3:  f  <  — “V^O+l)  -  p  <  1  .  (3-25) 
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Thus,  for  P  7^  3  Eq.  (3-22)  is  satisfied  in  one  range  of 
each  of  Eqs .  (3-16)  and  (3-17),  while  Eq.  (3-23)  is  satisfiea 
in  the  other  ranges.  The  regions  described  by  either  of 
Eqs.  (3-l6)  or  (3-1?)  are  consequently  possible  regions  of 
plastic  dissipation.  For  the  special  case  b  =  3  ,  (3-21) 

is  satisfied  in  the  range,  Eq.  (3-l8)*  branches  of  the 

curves  of  s^  and  t  against  U  ,  shown  in  Figs.  2  and  3, 
represent  therefore  valid  solutions. 


Equations  (3-l4)  and  (3-15)  already  give  s^  and  T 
as  functions  of  U  ,  while  the  remaining  quantities  J^,  u 
and  V  can  be  obtained  by  integration  from  Eqs.  (3-19^,l’,<l) , 
respectively.  To  determine  ,  Eq.  (3-19^)  1^-  inte¬ 

grated  by  parts 


J^(U) 

k 


(3-26) 


After  a  series  of  manipulations  shown  in  detail  in  Appendix 
A,  the  integral  in  Eq.  (3-26)  is  obtained  in  closed  form  in 
terms  of  elliptic  integrals, 


J  (u)  J  (U  ) 

- V—  = «(")  - 


(3-27) 


where  J  (U  )  is  the  value  of  J  at  the  boundary  U  of 
the  region.  When  0  <  U  <  V [3/3  the  function  H(u)  is 


+h(u)  = 


3.^ 

^  I  T 
V  I  X 


F(0,k)  +  V3(b+1)  E(0,k) 


+  ^  n  (0.F,k)  -  (b-3 


v^+1 


(3-28) 
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where  the  sign  of  h(u)  must  be  chosen  to  correspond  with 
that  of  s  in  Eq.  (3-14).  F,  E,  and  IT  are  the  elliptic 
integrals  of  the  firsts  second^  and  third  kind,  respectively, 
and 


sin  0 


SufT6+Tr~ 


(3-29) 


g  +4 
4(g+l) 


(3-30) 


When  V g/3  U <( n/  (4+g)/3  the  term  tanh“^  in  Eq.  (3-28) 
is  to  be  replaced  by  coth"^.  The  function  H(u)  is  well- 
behaved  except  for  U  =  \/g/3  where  the  tanh“^  (or  coth-^) 
term  becomes  infinite.  When  g  =  3^  the  term  does  not  appear. 

Similar  integrals  for  the  velocities  u  and  v  are 
also  evaluated  in  Appendix  A.  The  expression  for  u  in- 
volves  only  transcendental  functions,  while  the  one  for  v 
is  in  terms  of  elliptic  integrals,  except  for  g  =  3  ^  where 
the  elliptic  integrals  reduce  to  simpler  transcendental  func¬ 
tions  . 


B, 


DISSIPATIVE  SHOCK  FRONTS 


The  wave  propagation  problem  under  study  being  hyper¬ 
bolic,  the  possibility  of  discontinuities  in  the  solutions 
must  also  be  considered,  and  in  addition  to  dissipative  re¬ 
gions,  potential  locations  of  fronts  of  discontinuity  in¬ 
volving  energy  dissipation  must  be  explored.  The  general 
study,  [5],  of  elastic-plastic  wave  propagation  shows  that, 

in  an  xQcsax  c  xa  o  l.  xv^  "“pxao  j-v-  jl  jlu  a.  o  i-w  j '-x-.  w  uv-/  jj.v-.ju'-* 

condition  (2-1),  fronts  of  discontinuity  can  occur  only  in 
locations  where  the  shear  stress  tangential  to  the  front 
vanishes.  The  discontinuities  are  restricted  to  the  mean 
stress,  J^/3  j  and  to  the  particle  velocity  normal  to  the 
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front,  the  velocity  of  which  is  V K/p  ,  where  K  =  t3G/3  is 
the  (elastic)  bulk  modulus.  The  complete  solution  of  the 
problem  considered  here  being  a  function  of  U  as  defined 
t>y  Eq,  {3~^) )  such  a  front,  if  it  occurs,  is  located  at 


U  = 


(3-31) 


The  requirement  t  =  0  defin  s  the  value  of  s  at  the 
front 


s 

X 


k 


(3-32) 


while  the  discontinuities 
mentum  considerations 


AJ 

1 


and  Au  are  related  by  mo- 


AJ^  +  3PpG  Au  =  0  .  (3-33) 

In  addition,  the  condition  that  energy  at  the  front  is  dissi¬ 
pated  requires 

Ajys^  >  0  (3-3^) 

where  "A"  signifies  the  value  behind  the  front,  minus  the 
value  ahead. 

The  fact  that  a  discontinuity  in  the  solutions  may  oc¬ 
cur  at  U  V p/3  is  the  cause  of  the  breakdown  of  the  dif¬ 
ferential  relations  (3-19)  noted  previously. 

C .  NON-DISSIPATIVE  REGIONS  AND  NON-DISSIPATIVE  SHOCK  FRONTS 

In  non-dissipative  regions,  the  governing  equations  (2-15) 
and  (2-21),  when  written  in  terms  of  the  non-dimensional 
variables  (3-2),  and  of  the  independent  variable  U  become 
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while  the  yield  condition  Eq.  (2-22)  becomes 


(3-35) 


(3-36) 


Relations  (3-35)  being  a  set  of  homogeneous  linear  equations, 
it  follows  that 

s'  =  J‘  =  u'  =  t'  =  v'  =0  (3-37) 

^  1 


unless  the  determinant  of  the  system  vanishes,  i.e.. 


U(l-U2)  (4-i-p-3u2)  =  0  .  (3-38) 


Separating  each  non-dissipative  region  which  contains 
positive  roots*  of  Eq.  (3-38)  into  two  or  three  regions  ex¬ 
cluding  these  roots,  Eqs.  (3"37)  require  that  all  stresses 
and  velocities  within  the  regions  formed  in  this  manner  are 
constant . 


*  The  negative  values  of  U  are  of  no  consequence  because 
U  )>  0  ,  and  the  root  U  =  0  can  be  shown  to  be  trivial. 
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Without  further  analysis,  one  expects  in  a  non-dissipa- 
tive  region  the  discontinuities  known  from  the  theory  of 
plane  elastic  waves,  i.e.,  P-  and  S-fronts,  and  indeed  the 
positive  roots  of  Eq.  (3-38) 


U 


P 


(3-39) 


define  locations  moving  with  the  proper  velocities.  The 
solution  may  therefore  have  discontinuities  at  these  points. 
At  discontinuities  tv~  ,  s  ,  AJ  ,  and  Au  may 

occur,  related  by 


Aa  ;  As  :  AJ  :  Au  = 

X  X  1 

while  at  U  discontinuities 
S 

where 


^  *  4+p  •  •  '/4+F  (3  ^0) 

AT  and  Av  are  possible. 


AT  ^  -  N  j:G  Av  .  (3-^1) 

Obviously,  the  intensity  of  the  discontinuities  must  not 
violate  the  yield  condition  '3-36). 
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IV.  CONSTRUCTION  OF  SOLUTIONS 


In  the  previous  section,  solutions  of  the  basic  differ¬ 
ential  equations  have  been  obtained  for  regions  with  and 
without  energy  dissipation,  not  knowing  the  extent  of  these 
regions  in  the  complete  solution  to  be  found.  Further,  there 
are  three  potential  locations  of  discontinuities,  one  with, 
and  two  without  dissipation.  These  partial  solutions  must 
now  be  combined  to  form  a  complete  solution  satisfying  the 
externally  prescribed  conditions.  Due  to  the  multitude  of 
possible  combinations,  and  due  to  the  lack  of  a  general 
existence  and  uniqueness  theorem,  it  is  imperative  to  use  a 
systematic  procedure  in  the  construction  of  solutions  to  be 
certain  that  no  solution  is  overlooked. 


The  yield  condition,  Eq.  (2-1),  being  an  even  func¬ 
tion  of  the  components  of  stress,  the  complete  solutions  will 

also  be  even  in  a  and  in  t  .  Without  loss  of 

o  o 

generality,  the  following  will,  therefore,  consider  only  the 
case  of  compressive  surface  pressure,  <(  0  ,  and  positive 
shear,  0  <(  <(  k  ,  Consequently,  the  lower  sign  must  be 

taken  for  s  in  Eq.  (3-l^)>  and  correspondingly,  in  Eq. 
(3-28)  and  in  Appendix  A  for  the  functions  H(u)  ,  H^(u)  and 
H^(U)  . 

The  prescribed  surface  loading,  Eq.  (3-1)^  leads  to  the 
boundary  conditions  for  U  =  0  , 


Ox(0)  =  s^(0)  +  ^  J^(0)  = 
'r(O)  = 


=  o. 


(4-1) 
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where  the  nature  of  the  material  requires  !  ^ 

half-space  being  originally  stress  free  and  at  rest,  all 
stress  components  and  velocities  vanish  for  U  ->  oo  ,  furnish¬ 
ing  the  boundary  conditions 

lim  [s  ,J  ,u,v]  =  0  .  (^-2) 

Because  the  solutions  in  dissipative  regions  change  in 
character  depending  upon  the  value  of  p  ,  separate  treatment 
for  the  cases  f:'  )»  3  and  [3  \  3  required. 

A.  SOLUTIONS  FOR  3 

Starting  with  the  fact  that  for  U  oo  all  stresses 
vanish,  the  solutions  will  be  developed  by  consideration  of 
the  conditions  for  gradually  decreasing  values  of  U  .  The 
stresses  given  by  Eg.  (4-2)  do  not  satisfy  the  yield  condi¬ 
tion  (3-10),  and  the  first  changes,  from  vanishing  stresses 
to  initiation  of  yield,  must  therefore  be  elastic.  Such 
changes  can  only  be  discontinuities  at  U  =  Up  =  \/  (4+p)/3  and 

at  U  =  Ug  =  1  ,  described  by  Eqs.  (”^-40)  and  (3“4l)  . 

Case  A . 1  a  combination  of  discontinuous  changes 

Aa  and  At  in  their  respective  locations  are  selected 
so  that  yield  according  to  Eq.  (3“10)  is  not  reached,  no 
dissipative  regions  or  shock  can  occur  anywhere,  and  the 
stresses  introduced  at  the  P-  and  S-fronts  must  remain  con¬ 
stant  up  to  the  surface,  U  -  0.  The  resulting,  completely 
elastic  solution  leads  to  surface  loads 


If  the  combination  of  Aa  and  At  just  satisfies  the 

X 

yield  condition,  but  no  dissipative  region  is  utilized,  the 
above  solution  still  applies.  It  is  therefore  valid  for  the 
range 


0>^> 

^  -  2  v/I 


(4-4) 


Case  A. 2  Consider  next  the  situation  when  the  combina¬ 
tion  of  Ao  and  At  satisfies  the  yield  condition  (3-10), 
and  AT  0  .  For  ^  >  3  the  yield  condition  being  satis¬ 

fied  only  for  values  u  1  ,  the  plastic  shock  at 
U  =  V b/3  )>  1  cannot  occur,  and  of  the  two  potential  ranges 
for  dissipative  regions.  Eg.  (3-16),  only  one 

0  <  U  <  U  <1 
—  1  2  — 


remains  to  be  considered.  For  ^  =  3  >  the  value  of  the 
shear  t  =  At  7^  0  at  U  -  =  1  also  precludes  a  plas¬ 

tic  shock,  and  the  above  range  is  again  the  only  one  to  be 

considered.  At  the  terminal  point  U  of  the  dissipative 

2 

region,  if  any,  the  value  of  the  shear  stress  must  be  equal 
to  that  in  the  adjoining  non-dissipative  region,  t  =  AT  . 
Figure  3  shows  t  as  a  function  of  U  for  typical  values 
of  3  .  For  any  b  ,  there  is  one,  and  just  one,  value  u 

2 

which  corresponds  to  t  =  At  ,  where  a  region  of  plastic  dis¬ 
sipation  may  begin.  The  other  end  point,  U  <(  U  ,  of  the 

1  2 

dissipative  region,  may  be  selected  at  will,  so  that  it  may 

correspond  to  any  value  of  t^  >  At  .  As  approaches 

U  ,  the  value  t  approaches  At  ,  until  in  the  limit  the 
2  o 

solution  for  a  =  a  obtained  in  A,1  is  again  reached.  The 

o  1 

region  >  U  )>  0  is  necessarily  non-dissipative  with  con¬ 
stant  stresses. 
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The  entire  stress  history  is  obtained  by  computing  s 

from  Eq.  (3-l4),  t  from  (3-15),  and  J  from  (3-27)  and 

1  2^ 

(3-28).  The  normal  stress  at  the  surface.  o=s  (u  )+^  (u  ), 

is  a  function  of  ,  At  and  of  the  corresponding  values 

U  and  u  , 

1  2 


k 


(4-5) 


This  expression  for  depends  on  At  through  the  second 

and  fourth  terms  only,  both  contributions  having  the  same 
sign.  Both  terms  vary  monotonically,  and  because  in¬ 

creases  with  decreasing  At  ,  sj  1  -  (AT/k)^  and  ^(U^)  both 
vary  inversely  to  At  .  For  a  given  value  of  t^  ,  the 

largest  surface  pressure,  0  ,  is  therefore  obtained  for 

2 

At  0  .  Thus,  case  A. 2  applies  if 


-0  \  -0^  <  -a 

1  -L  o  A.  2 


(4-6) 


where  a  is  defined  by  Eg.  (4-4),  while  a  is  an  implicit 
1  2 

function  of  t_  and  6  , 
o 


— ^  ^  H(U  [t  ,p],p)  - 

The  velocity  history  may  be  obtained  from  Eqs . 
and  (3-41),  and  from  Appendix  A. 


J  H(l,6)  . 

(4-7) 

(3-40) 


Case  A.^  Finally,  let  the  discontinuity  Ac.,  be  '^uch 
that  the  yield  condition  will  be  satisfied  by  the  discontin¬ 
uity  Ac  at  U  =  Up  alone, 

X  ir 


Ao 

“k 


X  ^  _  4  +  B 
2  sf3 


(4-8) 
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both  of  the 


The  yield  condition  being  satisfied  at  U  ^  Up  ^ 
two  dissipative  regions,  Eq.  (3-16),  must  be  considered. 
However,  the  shear  t  just  behind  the  P-front  being  zero. 
Fig.  3  shows  that  no  dissipative  region  can  actually  occur 
for  U  ^  n/pTB  because  the  starting  and  end  points  of  such 
a  region  would  be  identical,  =  U  =  V p/3  .  The  stresses 
in  the  range  Up  )>  U  )>  \^3/3  can  therefore  not  vary  and  are 

given  by  Eq.  (4-8).  At  U  =  Vp/3  ,  the  necessary  condition 
T  =  0  being  satisfied,  a  dissipative  shock  may  occur.  The 
strength  of  this  shock,  described  by  AJ^  ,  is  arbitrary  ex¬ 
cept  for  the  sign, 

AJ^  <  0  .  (4-9) 

In  the  range  V P1/3  )>  U  )>  1  ,  there  is  no  alternative  to  a 

non-dissipative  region  of  constant  stress.  At  U  =  1,  the 

yield  condition  excludes  the  possibility  of  an  elastic  change 

in  shear,  but  Fig.  3  indicates  that  a  dissipative  region  is 

possible  from  U  =1,  corresponding  to  x  =  0,  to  any  de- 

2 

sired  value  U^  ,  corresponding  to  the  surface  shear  x^  . 

The  surface  pressure  is  given  by 

O  =0  +  AJ  .  (4-10) 

o  201 

If  AJ  is  selected  equal  to  zero,  the  solution  obtained 
1 

is  identical  to  that  obtained  in  Case  A. 2  for  a  =  a  .  By 

o  2 

using  sufficiently  large  values  of  j  AJ  |  ,  one  obtains 

solutions  for  any  value  (a  I  )>  I  a  I  of  the  surface  pres- 

O  2 

sure  without  upper  bound.  Equation  (4-10)  defining  AJ^ 
uniquely  as  a  function  of  I  I  )>  I  a  I  ,  the  solutions  in 
the  range  are  again  unique,  and  do  not  overlap  Case  A. 2. 
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The  corresponding  velocity  history  may  be  obtained  from 
Eq.  (3“33)  in  addition  to  the  relations  noted  in  A. 2. 

The  three  cases  considered  exhaust  all  the  possible  com¬ 
binations  of  regions  and  fronts  for  3  >  3-  There  is  one  and 
only  one  solution  for  any  possible  combination  of  surface 
loads  and  .  The  regions  in  the  ,  T^-plane 

governing  each  case  and  the  separation  lines  and 

are  shown  in  Fig.  4^  drawn  for  P  =  5  •  The  appropriate  ex¬ 
pressions  for  the  solution  are  listed  in  Tables  I  to  III. 
Typical  distributions  of  and  t  for  the  non-elastic 

cases  are  shown  in  Figs.  5  ^rid  6. 

B .  SOLUTIONS  FOR  b  <  3 

Reasoning  exactly  as  for  p  )  3  ^  the  stresses  must 
vanish  for  large  values  of  U  until  the  P-  or  S-fronts  are 
reached,  and  elastic  discontinuities  in  stress,  Ao  and  At, 
occur  to  initiate  yield. 

Case  Bol  If  the  discontinuities  Aa,.  and  Ax  are 
sufficiently  small  so  that  yield  is  not  reached,  the  situa¬ 
tion  is  exactly  as  in  Case  A.l  for  |3  /  3  ^nd  Eqs .  (4-3) 
and  (4-4)  apply. 

Case  B.2  If  the  stress  discontinuities  Aa  and  Ax 
-  x 

are  such  that  the  yield  condition  is  just  satisfied,  but 

Ax  0  ,  a  dissipative  region  or  shock  is  impossible  unless 

U  <(  1  .  Due  to  Ax  0  ,  the  dissipative  shock  at 

U  =  sfJ73  <(  1  cannot  occur,  since  it  requires  t  =  0  ,  but 

a  dissipative  region  may  start  at  a  point  U  <(  V p/3  corre- 

2 

sponding  to  x  =  Ax  ,  Continuing  this  region  to  a  point 
U  <)  U  corresponding  to  x  =  x  >  Ax  ,  a  solution  is  ob- 
tamed  for  a  value  ,  not  known  a  priori.  As  in  Case  A. 2, 

the  entire  stress  history  is  obtained  from  Eqs.  (3-14),  (3-I5), 
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(3“27)j  and  (3-28),  while  the  normal  stress  at  the  surface, 

a^,  is  given  by  Eq.  (4-5).  The  dependence  of  on  At 

is  again  monotonic  so  that  the  solutions  found  in  the  range 

are  unique.  It  is  important  to  note  that  the  function  H(u) 

becomes  infinite  as  U  approaches  V p/S  .  Therefore,  for  a 

given  ,  it  is  evident  that  in  the  limit  At  ^  0  , 

U  Vp/3  ,  and  the  surf" ce  pressure  will  increase  without 

bound,  a  ->  -oo  . 
o 

The  solution  described  applies,  therefore,  for  any  value 
of  the  surface  pressure  not  covered  by  Case  B.l,  except  that 
the  condicion  }  Ar  ^  0  excludes  the  case  =  0- 

Case  B . ^  The  remaining  possibility  is  a  discontinuity 
Aa  at  U  ,  sufficient  to  satisfy  the  yield  condition.  In 

^  Jr* 

the  adjoining  non-dissipative  region,  a  =  Ao  and  t=At=0, 
which  permits  a  dissipative  shock  of  arbitrary  strength 
AJ  <(  0  at  U  =  V p/3  .  This  shock  may  be  followed  by  a  non- 
dissipative  region  of  constant  stress,  furnishing  solutions 
for  T^  =  0  for  any  value  of  a  not  covered  by  Case  B.l. 

If  one  attempts,  starting  as  in  the  previous  paragraph, 

to  employ  a  dissipative  region  behind  the  shock,  i.e.,  from 

U  =  V 5/ 3  to  some  value  U  <(  U  ,  the  singularity  in  h(u) 

2  1  2 

at  V p/3  leads  to  =  -cx)  ,  a  physically  unrealistic  situ¬ 

ation  . 

The  three  cases  considered  furnish,  therefore,  one  and 

only  one  solution  for  each  combination  of  o  and  t 

o  o 

The  regions  in  the  ,  T^-plane  governing  each  case  are 

shown  in  Fig.  7-  The  appropriate  expressions  for  the  solu¬ 
tion  are  listed  in  Tables  I,  II,  and  IV.  Typical  distribu¬ 
tions  of  0^  and  t  for  Case  B.2  are  shown  in  Fig.  8.  The 
velocity  history  for  (:'  <(  3  n\ay  be  obtained  in  a  manner 
analogous  to  that  for  (3  >  3- 
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V.  SUMMARY  OF  RESULTS  AND  DISCUSSION 


The  stresses  and  velocities  in  an  ideal  elastic-plastic 
half-space,  subject  to  the  v.  Mises  yield  condition,  have 
been  obtained  due  to  normal  and  shear  stresses  on  the  sur¬ 
face,  both  of  which  are  step  functions  of  the  time.  The 
solutions  are  in  closed  form,  containing  elliptic  integrals 
and  simpler  transcendental  functions.  Expressions  for  the 
stresses,  including  pertinent  basic  relations,  are  tabulated 
in  Tables  I  to  IV,  while  the  expressions  for  the  velocities 
are  given  in  Appendix  A.  The  results  indicate  a  number  of 
ranges  where  different  solutions  apply,  depending  on  the  in¬ 
tensity  of  the  surface  loads,  and  on  the  value  of  a  material 
parameter  f  '^^ich  varies  with  Poisson's  ratio. 

The  dependence  of  these  ranges  upon  the  applied  surface 
stresses  o  and  t  is  shown  in  Fig.  4  for  B  z*  3  ^  and 
in  Fig.  7  for  B  s  3 • 

Excluding  the  cases  where  and  are  so  small 

that  no  plastic  deformation  occurs.  Figs.  6,  and  8  show 

the  stresses  for  a  typical  case  of  each  of  the  m^^jor  situa¬ 
tions.  Figures  5  and  6  apply  if  B  .  3  ^  which  corresponds 

to  Poisson's  ratio  v  ^  1/8.  The  non-dimensional  ordinate 
U  in  the  plots  is  proportional  to  x/t  so  that  the  fig¬ 
ures  illustrate  the  x-distribution  of  the  stresses  at  any 
Lime  L.  For  an  applied  compressive  stress,  cIq  \  0  ,  rn  the 
range  ^  ^  ”*^0  ^  ”"^2 ^ ^  ^  stress  plot  is  given 

by  Fig.  5-  There  is  an  elastic  precursor  consisting  of  a  P- 
wave  of  uniform  strength,  followed  by  an  S-wave,  also  of 
uniform  strength,  so  that  the  combination  of  the  two  states 
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of  stresses  just  satisfies  the  yield  condition.  Proceeding 

towards  the  surface  there  is  no  plastic  deformation  or 

energy  dissipation  until  a  point  U  =  U  is  reached.  Plastic 

2 

deformation  accoi  panied  by  energy  dissipation  occurs  in  a 
region  from  U  to  .  In  this  region,  and  t  increase 

until  the  surface  values  a  and  r  are  reached.  Finally, 
betwe*' a  the  end  of  the  plastic  region  at  U^,  and  the  sur¬ 
face,  the  stresses  are  constant.  If  -a^  )>  -o,  ,  Fig.  6,  there  is 
again  an  elastic  precursor  consisting  of  a  P-wave  of  uniform 
strength  so  that  the  yield  condition  is  just  satisfied. 

At  U  =  V p/3  there  is  a  plastic  shock,  a  discontinuity  in 
a  without  change  in  r  ,  followed  by  a  region  of  constant 
stress.  From  U  =  1  to  U  =  U  ,  there  is  a  plastic  region 

where  a  and  t  increase,  just  as  in  Fig.  until  the 
yc 

surface  values  o  ,  t  are  reached.  For  0  <(  U  <(  U  , 

o  o  1 

there  is  again  a  region  of  constant  stress. 

For  P  <  3  ,  V  <  1/8  ,  Fig.  8  applies  for  -o^>-a^(T^,p) 

and  '’’q  ^  O’  stress  history  in  Fig.  8  is  similar  to 

Fig.  5  3,3^  there  is  an  important  difference. 

As  the  end  point  U  of  the  plastic  region  approaches 
_  2 

vp/3  j  che  surface  pressure  -o^  increases  beyond  any  bound, 
so  that  a  solution  with  a  plastic  shock  never  occurs,  ex¬ 
cept  if  T  =  0  . 

^  o 

In  view  of  the  fact  that  no  general  existence  and 
uniqueness  theorem  for  the  problem  of  wave  propagation  in 
elastic-plastic  media  is  available,  it  i  pertinent  to  note 
that  solutions  v/ere  found,  for  any  combination  of  Poisson's 
ratio  V  ,  surface  pressure  “Oq  ’  surrace  shear  I  ^  ^  ^ 

and  that  these  solutions  are  unique. 

While  the  method  of  solution  used  in  Sec.  Ill  does  not 
appear  to  utilize  the  characteristics  of  the  system  of 
partial  differential  equations,  the  characteristic  velocities 
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appear  in  the  analysis  in  a  disguised  form.  Inverting  Eqs . 

15)1  so  that  the  non-dimensional  variable  U  , 

Eq.  (3~'4),  becomes  a  double  valued  function  of  the  state  of 
stress,  it  IS  shown  in  Appendix  B,  that  the  characteristic 
velocities  are  related  to  the  values  of  U  , 


Figures  2  and  3  demonstrate,  therefore,  the  dependence  of 

the  characteristic  velocities  V.  on  the  values  of  s  or  t. 

1  X 

The  availability  of  the  characteristic  velocities  , 

and  of  the  relations  along  tiie  characteristics,  also  given 
in  Appendix  B,  are  an  invitation  to  their  use  for  numerical 
analysis  in  the  case  of  general  surface  stresses  and 

T^(t).  However,  the  moving  boundaries  between  dissipative 
and  non-dissipative  regions  pose  difficulties  requiring 
further  study  as  outlined  in  Appendix  B. 

The  approach  used  in  Sec.  Ill  is  readily  applicable  to 

other  types  of  yield  conditions.  The  case  for  the  Tresca 

condition  becomes  slightly  more  involved  than  for  the  v. 

Mises  condition,  since  there  are  several  sub-cases  depending 

on  whether  or  not  is  one  of  the  extreme  principal 

stresses  entering  the  yield  condition,  I  c  -  0  .  I  =  2k. 

^  ^  max  min 

In  the  sub-case  which  controls  in  the  present  problem,  o 

z 

is  th'T'  intermediate  principal  stress  and  the  relations  be¬ 
tween  s  or  i  and  U  ,  Eqs.  ^J-ig,  15) j  change,  but 
without  affecting  the  nature  of  the  solutions  and  with  only 
minute  numerical  effects. 
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Another  potential  generalization  of  the  present  paper 
is  a  two-dimensional  problem,  the  determination  of  the  steady- 
state  solution  for  an  elastic-plastic  half-space  subjected 
to  a  step  pressure  on  the  surface  x  =  0  ,  moving  with  super- 
seismic  velocity  .  Using  a  system  of  polar  coordinates 

moving  with  the  front  of  the  load,  Fig.  9,  dimensional  con¬ 
siderations  require  the  stresses  to  be  functions  of  the 
angle  9  ,  but  not  of  the  radius  r  .  This  permits  the  re¬ 
duction  of  the  partial  differential  equations  into  a  set  of 
ordinary  ones,  containing,  however,  a  larger  number  of  un¬ 
known  dependent  variables  than  in  the  problem  treated  here. 
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APPENDIX  A 


INTEGRATIONS  REQUIRED  FOR  THE  DETERMINATION 
OF  J  ,  u ,  and  v 


1.  EVALUATION  OF 

Starting  with  Eg.  (3-26)  in  the  text,  substitution  of 
the  value  of  s  /k  from  Eq.  (3-14)  yields 


^  _  2ue  v/T  ./4+6-3U^ 

k  3U^  -  P  ]/  ^  + 


.  f  .  /u'^(4+a-3U^)  2UdU  , 

+  66  nA3  J  1/  (3U^-r-^  1  const, 


(A-1) 


requiring  the  evaluation  of  the  integral 


U 


I(U)  = 


/HS+OIir  2fde 


u® 


6  + 


_ , _ ^  =  f  ./n(4-f6-3Ti)  .  _dr] _ ^ 

3^^-3)‘^  J  y  6  +  r,  (3t]-6)-^  • 

o 

(A-2) 


Multiplying  numerator  and  denominator  in  the  last  integral 
by 


Vri(4+r-3Ti) 


l(u)  becomes 


-  n) 

5^2 


I(U)  = 


dT] 


i  /)4^- 


(A-3) 


t))t](t)+6) 
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which  is  of  the  form 


f 


r( t)dt 


>/  ( a-t )  ( t-b )  ( t-c ) 


where  R(t)  is  a  rational  function  of  t  .  This  integral 
may  be  converted^  Byrd  and  Friedman  [6],  No.  235* 20^  with 
a>y>b>c  ,  to 


/■ 


r(  t)dt 


V ( a-t) ( t-b ) ( t-c ) 


=  g 


where* 


sn^u 


a-c ) ( t-b 


a-b) ( t-c 


P  = 


g  = 


a-b 

a-c 


V  a-< 


0  =  am  u 


snu  = 


1 

sin  0 


•  -1  1  a-c ) ( y-b ) 
V(a-b)(^:c) 


(A-5) 


and 


a  =  ,  b  =  0  ,  c  =  -(3  and  y 


=  U" 


(A-6) 


while 


r(t))  = 


,  2i4-R 


3  V  3  (t)  -  ^) 


^\2 


*  Note  the  difference  in  meaning  of  the  symbols  U  and  u. 
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Substituting  these  values  into  Eq.  (A-5) ,  one  has 


sn^u 

g 


4  ( b-H  )  ri 

(3+4) (b+n) 

:5+4 

"  4(3+1) 


sin  C 


snu 

1 


,  /  4(;-  +  l)U^~ 


(A-7) 


The  argument  of  R  which  appears  on  the  right  side  of 
Eq.  (A-4)  is 

3(  !H4  )  sn^u 
4(3+1)  dn^u 

and  the  integral  (A-3)  finally  becomes 


u 


I(U)  =  C3+4)^  r  ^  sn^u[4(3+l)dn^u  -  33  sn^u] 

'  T  J  [4(3+l)dn2u  -  3(3+4)sn2u]2 


3  ^3^ 


(a-8; 


Noting  the  relations  (Pierce  [7]:.  Nos.  730^  731) 


sn^u  +  cn^u 


=  1 


dn“u  -t-  k^sn-u  =  1 


Eq.  lA-8)  reduces  further  to 

u , 


I(U)  - 


lii±4I 


^  2 


43(3+1)' 


sn^u  cn^u 


(1  -  sn==u)- 


du 


(A-9) 
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Applying  [6J  No.  362. l8  for 


_  fi±i  v2  ^  3+4  ^  qS 

3+1  ^  4( 3+1)  4 

gives 


I(U)  =  n(0,a2,k)  -  E(0,k)  -  F(0,k) 


+ 


snu  cnu  dnu 


(A-10) 


where 


u 

E(u) 
n  ( u  j  ) 
and 


F(0,k)  =  elliptic  integral  of  the  first  kind 

E(0jk)  =  elliptic  integral  of  the  second  kind 

n( 0ja^ ,k)  =  elliptic  iricegral  of  the  third  kind 


snu  =  sin  0  = 
1 


4(3+1) 

(3+4)(3+U=") 


cnu  =  cos 
1 


^  _  -\/pL4±3;3u.^1-  ) 

“  V(3+4)(3+u^) 


dnu  =  J 1  -  k‘=^sin‘^0  ^  V 

1  p+u 


2  •  2. 


-  V— ^ 

V  0_LTT 


J 


(A-11) 


Using  Eg.  (A-ll),  the  last  term  in  Eg.  (A-10)  becomes 


snu  cnu  dnu 
- i_a. - i. 

^_2,. 


3+1^ 


23U 


n  tv2 


4  +  3  -  iiF 

i,p-rx;  VP^u-  ; 


(A-12) 


Returning  to  the  original  expression  (A-1) ,  noting  that 
1(0)  =  0  and  that  expression  (A-12)  cancels  with  the  first 
term  on  the  right  side  of  Eg.  (A-1)^  one  has 
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=  H(U)  +  const. 


(A-13) 


where 


k 


H{U)  =  +  n/3('p+i)" 


^  F(0,k)  +  E(0,k)  -  |('|^jn(0.a^,k)! 


J 


.  (A-14) 


If  the  constant  of  integration  in  Eq.  (A-I3)  is  expressed 

by  the  value  of  J  /k  at  the  end  of  the  region  U  ,  then 

1  2 


J  (U)  J  (U  ) 

- =  H(U)  -  H(U^)  . 


(A-15j 


For  any  given  value  of  U  ,  in  a  dissipative  region,  0 
may  be  found  from  relation  (A-ll)  while  the  F  and  E  func¬ 
tions  may  be  taken  from  any  of  several  tables  of  ell-ptic 
integrals,  e.g.,  [8j.  The  integral  of  the  third  kind,  II  , 
is  tabulated  in  Selfridge  and  Maxfield  [9]  for  various  values 
of  k  but  only  for  <(  1.  "Additional  formula"  on  page  xi 
relates  the  tabulated  values  with  those  for  )>  1  required 
here , 


n(<^,a",k)  +  n(c,^,k)  =  F(o,k) 


a 


(1-aM 


tan 


-1 


Since  3  1  and  >  k^  ,  each  of  the  radicals  above  will 

be  imaginary.  Equations  6  and  7  P^ge  51  i"  [10]  may  be 
used  to  obtain  the  last  term  of  Eq.  (A-16)  as  a  real  function 
of  a  real  argument. 
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Let  the  value  sought  be 


f  =  j  tan"^  iX  (A-'17) 

where  x  is  real  and  either  ^  1.  Inverting, 


y  =  i  tan  i  iIj  -  —  s in ( x+  i y )  _  sinh  v  cosh  y-i  sin  x  cos  x 
"  i  ^  ”  i  cos(x+iy)  "*  cos^x  cosh^  y+sin^x  sinh^  y 

(A-18) 

where  i^  =  x  +  iy.  The  condition  that  X  be  real  implies 
that  either 


7T 


X  =  n?:  or  x  =  (2n-l)  ^ 

which  in  turn  gives 


(A-19) 


f  =  tanh“^X  -  inTT  or  ^  =  coth”^X  -  i(  2n-l)^.  (A-20) 

In  either  case  the  imaginary  constant  may  be  included  as  part 
of  the  constant  of  integration  so  that  H(u)  remains  real. 

Using  the  appropriate  values  of  0,  and  k  ,  the 

formula  (a-i6)  for  11(0, a^,k)  becomes  for  U  <(  Vp73 


n(0,a^,k)  =  F(0,k)  -  n(0,^,k)  +  — - — V P+1  tanh  ^ 

^  V3 


3u  b+_ui_ 
^  V4+3-3u^ 


{A-21) 

and  Eq.  (A-14),  completely  in  terms  of  tabulated  functions,  is 
■fH(u)  =  F(0,k)  +  V3W1T  E(0,k)  +  3...^ip-3)  n(0,|-,k) 


V  p+i 


2 


4  VpT+T 


1 

M  _ 

/3  + 

mm 

y  4+p--!U‘^ 

(A-22) 
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where 


sin  ■h> 


(A-^,3) 


xliL 


4 1  1  /' 


For  U  V' iV3  ,  Eqs.  (A-21)  and  ';a-22)  must  be  modified  by 
replaciny  tanh-^  by  coth~^  .  The  function  H(U)  is  plotted 
in  Fig.  10  for  r  =  5  ^  3  iri  the  range  0  <(  U  ^  1  ^  and  for 
=-  2  <(  3  in  the  range  0  <(  U  \  s''  2/3  .  No  plot  for 
U  ^  V t/3  is  shown  because  this  case  does  not  actually  occur 
in  the  solutions  constructed  in  Sec.  IV.  An  infinity  arises 
at  U  =  '32/3  in  Fig.  10,  for  the  case  r  =  2  <  3,  because 
the  argument  of  tanh”*  in  F.q.  (A-22)  approaches  unity  and 
the  function  becomes  infinite.  In  the  case  P  3  3  ^  simi¬ 

lar  situation  exists,  but  on  the  branch  U  ,,  V p/3  which  is 
not  plotted. 

For  the  special  case  t  =  3  >  Eq.  (A-22)  becomes 


+H(U)  =  2  2  V3  E  ^3,-^)  (A-24) 

which  is  al'ways  finite. 

Figure  10  shows  typical  plots  of  H(U)  in  the  region 
where  Eq.  (A~22)  applies.  For  a  considerable  part  of  the 
range  H(U)  is  nearly  linear  ,  so  that  the  first  term  of 
the  power  series  for  H(Ui 

+  H(U)  -  2U  -]/^  (p+4)  (A-25) 

may  be  used  as  an  approximation. 


-39- 


2. 


EVALUATION  OF  u 


Differentiation  of  Eq.  (3-1^)  yields 


s  • 

X 

k 


6(4+6)  -  66U^  -  3U^ 

o  +V)^ 


(A-26) 


which,  together  with  Eq.  (3~19b),  gives  the  derivative  of  u 


/ 


±y 


3(e+u 

4+I3-3U 


6(4+B)  -  6BU^  -  3U* 

(P-SU'^XB+U^)” 


(A-27) 


and  by  integration 


+ 


u(u)  =  H^(U)  t  const. 


(a-28) 


where 


H  (U) 


6(4+6)  -  66U^  -  3U^ 

(mu^TTp+u^P^ 


2UdU.(A-29) 


Introduction  of  the  new  variable 


{A-30) 


gives 


When 


H  (U) 
u  ’ 

U  <  'JWl  .  i 

H^(U)  =  2 


36  +  66.3^  -  (4+6)2^^ 
(6-z^ ) ( 1+z^ ) 


dz  . 


.e 


z  <(  i  the  integral  becomes 


tan“^z  -  tanh 

2  -fi  J1 


(A-31) 


(A-32) 
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Reintroducing  the  substitution  (A-30)  gives 


2  tan--' 


2  Vt 


p(4+r^-3u^ 


_  ^^3  ■  M+t'.~3y^ 
2  V  i^-t 


(A-33) 


Expressing  the  constant  in  Eq.  (A-28)  by  the  value  of  u  at 

U  yields  finally 
2 


u(U) 


H  (u  )  (A-34) 
^  2 


where  the  sign  of  u  corresponds  to  that  of  s^  . 

It  should  be  noted  that  the  argument  of  tanh"'^  approaches 
unity  as  U  V p/ 3  .  Thus  behaves  similar  to  H(u) 

and  becomes  infinite  at  U  =  •'fW/3  >  unless  3  =  3- 

3 .  EVALUATION  OF  v 

Equation  (3~19d)  gives  v'  in  terms  of  s  ,  and  s'  . 

yi  ^ 

Combining  Eqs .  (3-14),  (3-15)-  ^nd  (A-26),  one  obtains 


/  c  +  U- 
V  ( i-u'')  (^-30""  ; 


b(44-t)  -  6bU^  -  30^^! 


(A-35) 


and  by  integration 


"I' 


y  A. 


=  H  (U)  +  const, 
v 


I  a  --5  A  \ 

V  / 


where 


U 


H  (U) 

V  ^ 


I 


(l-^‘^)(t-3^^ 


-.6^ 


e  -  36- 


.(a-37, 
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i 


Using  the  value  of  v  at  the  boundary  of  the  region  to 
evaluate  the  constant ^ 

±')l^  v(u)  =  ±‘]/^  v(u^)  +  H^(U)  -  H^(U^)  {A-38) 

where  the  sign  of  v  is  chosen  to  correspond  to  that  of  t  . 
Making  the  substitution  ^  in  the  integral  gives 


f  \l 

"p(4+i3)  -  o^Ti  -  3r|2 

J  V  (i-n)  (p-3t])t^ 

_  WW  J 

dri 


4+3  )  -  6^r]  -3i1 

/ 


2 


o  J  (I  -  n)  (1-1^1  )Ti(n+B) 


dTi  . 


(A-39) 


The  la‘‘t  expression  may  be  compared  with  No.  25^.41  in  [6], 

y 


u 


/ 


R(tidt 


V ( a-t ) (b-t ) ( t-c ) ( t-d) 


=  g 


r\  ^ 


c-g^d  srPu 
lii^sn^u J 


du  (a-40) 


where  R(t)  is  any  rational  function^  where 


and 


a>b>y>c>d 


sn^u 


(b-d)( t-c) 
(b-c)(t-d) 

(p.sU^rAl 

(a-c) (b-d) 


(A-41) 


(A-42a) 

(A-42b) 
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a-c ) (b-d ) 


(A-42c 


0  <  <  k- 


(A-itad 


sin  <t> 


(b-d) ( v-c^ 
(b-c j (y-d 


(A-42e ) 


The  cases  3  ?  ^  3  j  p  =  3  niust  be  treated 

separately. 

Case  3. a  When  (  3  j 

a=|,  b=l,  c=0,  d=-p,  y=u2  (A-43) 

satisfies  the  condition  (A-4lj  for  0  \  U  <(  1  ,  and  Eqs . 
(A-42)  give 


= 


p+l  ^ 


.  §  = 


2./3_. 

pO+1) 


(a-44: 


sin  <!)  =  snu  =  U 
1 


'g_±l 

p+lF 


(a-45) 


The  argunent  of  R  which  appears  on  the  right  side  of 
Eg,  (A-40j  is 


p 

sn^u 


(  1  -  77;^  sn^u, 


and  the  integral  (A-39)  becomes 


H  (U 

V 


V  P  (  pd- 1  ) 


8sn^u  +  —77-  sn'^u 

_ 1  1 _ 

1  p  ^ 

■  "TTf  sn'^u 


du.  i'A-46 
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The  use  of  the  identity  dn^u  =  1  -  k^sn^u  ,  and 
337.01,  3^'2.11  in  [6],  together  with  Eq.  (A-44), 
the  final  form  of  the  integral 


Nos.  336.01, 
results  in 


H„(u)  =  37^  -  3,y^  ^(0,s^k) 

(a-47) 

/\  /\  yS 

where  k,  a,  and  0  are  given  by  Eqs .  (A-44)  and  (A-45) . 
(in  contrast  to  the  situation  encountered  when  determining 
H(u),  being  here  less  than  unity,  the  function  11  is 

tabulated . ) 

Case  3.h  When  3  3^  the  inequality  (A-41)  is  satis¬ 

fied  in  the  range  0  <(  U  <(  s/3/3  t)y 


a  =  l,  b=|-,  c  =  0,  d  =  -3,  y  =  u^ 
Eqs.  (A-42)  give 


and 


(A-48) 


(A-49) 


The  argument  of 


sin  0  =  sn  u  =  -  .  (A-50) 

R  on  the  right -hand-side  of  Eq.  (A-40)  is 


£ 

4 


(1 


2 

sn^u 


and  the  integral  (A-39)  becomes 
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1 


u . 


H  (U)  =  — ^ 

^  2/F 


/ 


4-rt  -  (2tr'4l)sn^u  -  sn^u  dn^u 


1  2  ^ 
1  -  jj-  sn^L 


du  .  (A-51) 


Using  Nos.  336.01,  337-01,  and  362.11  in  [6]  gives 


H^(U)  -  2  Ft  0,k)  +  — ^  E^0,k)  -  ^  v'“  n(  ,  ,^,k)  {A-52) 

/s 

where  k  and  0  are  given  by  Eqs .  (A-49)  and  (A-50),  re¬ 
spectively  . 

Case  3.C  When  r  3  -  Eg.  (A-39)  reduces  to 


Hv(U) 


(7+n)dq 
(3+’i)  t(3+3J 


(A-53) 


which  may  be  evaluated  by  using  the  substitution 


The  final  result  is 

H^(U)  =  tanh-^  -  -  IA-54) 

v’  3+u  \  3  3‘^0'“ 

The  same  result  could  also  be  obtained  by  substituting 
^  =  3  into  Eqs.  ,A-44)  and  .A-47)-  or  Eqs.  (A-49)  and 
(A-52).  For  the  resulting  value  of  the  modulus,  k  =  1  ,. 
the  elliptic  integrals  reduce  to  simpler  functions. 


APPENDIX  B 


USE  OF  THE  METHOD  OF  CHARACTERISTICS  WHEN  THE 
APPLIED  SURFACE  STRESSES  ARE  TIME  DEPENDENT 


In  Sec.  II,  the  basic  partial  differential  equations 
governing  propagation  of  plane  waves  in  both  dissipative 
and  non-dissipative  regions  were  obtained.  In  this  Appendix 
it  will  be  demonstrated  that  the  characteristic  velocities, 

,  of  these  differential  equations  may  be  obtained  from 
Eqs.  (3-l4,  15)  relating  s  ,  r  and  U  ,  but  that  diffi- 
culties  arise  in  the  application  of  the  method  of  charac¬ 
teristics. 

1.  DISSIPATIVE  REGIONS 


Equations  (2-15),  ( 2-1? )> and  Eq.  (2-19)  differentiated 
with  respect  to  time,  using  also  Eq .  (3“3)^  give  the  quasi- 
linear  set  of  differential  equations 
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Using  matrix  notation,  Eg.  (B-1)  becomes 


A'u  +  b'u^  -  0  (B-2) 

X  t 

where  u  is  the  column  matrix  of  the  unknowns,  and  I  b'I  i  0, 
Premultiplying  Eg,  {B-2)  by  b'"^  gives 


Au  +  u. 

X  t 


0 


(B-3) 


where 


A  =  b'“^  a' 


(B-4) 


Using  the  procedure  in  [11],  the  characteristic  velocities 
V  are  obtained  from 


A  -  VI  '  =  0 


(B-5) 


where  I  is  the  unit  matrix,  and 
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Making  the  substitution 


V 


U 


the  condition  (B-5)  gives 


u2[(i-u2)0-3u2)  -  (^)^0+u2)]  =  0  . 


(B-7) 


(B-8) 


This  equation  is  identically  satisfied  by  Eg.  (3“15)  so 
that  tl j  values  U  as  functions  of  t  (or  s  )  define  the 
non-vanishing  characteristic  velocities. 

Solution  of  Eq.  (B-8)  for  gives,  in  addition  to 

the  root  =  0  ,  always  two  different  positive  roots 

so  chat  there  are  four  characteristic  velocities,  V  =-V  , 

1  a 

V  =  -V  ,  in  addition  to  the  degenerate  double  value 

3  4 

Vo  =  0. 

After  appropriate  manipulation,  the  compatibility  equa¬ 
tions  along  the  four  non-vanishing  characteristics  become 


As  +  ^  AJ 
X  3  1 


pV^Au  + 
pVj^^  - 


k2 


S  T 
X 


3k^ 


3^ 


^  _ 
3 


Av  =  0 


(B-9) 


where  As  ,  AJ  ,  etc.,  are  the  increments  of  the  respective 
functions  along  the  characteristics.  The  degenerate  double 
value,  =  0^  furnishes  the  relations 


4  ®x^®x  +  =  0  (B-10) 

2  S  GA  =  AJ  -  AS^  .  (B-11) 

X  3b  1  X  ^  ^ 
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The  first  of  these  equations  can  bo  recognized  as  the  in¬ 
crement  of  the  yield  condition. 

The  five  Eqs .  (B-9.  10;  permit  the  step  by  step  evalu¬ 
ation  of  the  five  quantities  s  .  J  ,  u,  and  v  within 

X  1 

a  dissipative  region,  while  Eq.  .B-llj  provides  the  neces¬ 
sary  check  that  a  >  0,  required  in  a  dissipative  region. 

2.  NQN-DISSIPATIVE  REGIONS 

Starting  with  Eqs.  (2-15)  ^nd  (2-21),  one  finds  the 
equation 

V[V"  -  -  ^]  =  0  (B-12j 

P  o  P 

for  the  characteristic  velocities.  There  is  one  degenerate 
value 


and  the  four  well-known  values 
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2,3 

- 

At  + 
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(B-17) 
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These  relations  are  applicable^  provided  the  increments 
of  the  stresses  are  such  that  the  inequality  (2-22)  remains 
satisfied . 

The  appearance  of  the  degenerate  value  Vq  =  0  is  due 
to  the  fact  that  the  formulation ^  Eqs.  (2-15)  and  (2-2l)j 
uses  s^  and  J  as  separate  variables ^  instead  of  the 
single  unknown  a  =  s  +  J  /3  sufficient  for  the  solution 
of  a  purely  elastic  problem. 

3 .  DISCONTINUITIES 

Elastic  discontinuities  of  well  known  nature  may  propa¬ 
gate  with  the  velocity  of  P-  or  S-waves,  +  Cpj  +  Cg  ,  re- 
spectivelyj  while  dissipative  shock  fronts,  discussed  in 
Subsec.  III.B,  have  a  velocity  ^  =  +  Vk/p  ,  and  arv;  restricted 
to  locations  where  t  =  0  . 

When  constructing  solutions  involving  discontinuities, 
analytical  relations  are  required  to  follow  the  interaction 
of  elastic  and  plastic  fronts  with  adjoining  continuous 
plastic,  and  elastic  regions,  respectively.  There  is  no 
particular  difficulty  in  deriving  the  appropriate  relations. 
They  were  not  derived  here,  since  they  are  not  required  for 
the  discussion  which  follows . 

if .  DIFFICULTIES  IN  THE  DETERMINATION  OF  SOLUTIONS 

The  difficulties  encountered  are  best  illustrated  by 
considering  the  Cauchy  boundary  value  problem.  The  usual 
finite  difference  technique  can  be  used  in  any  location 
which  is  known  to  be  entirely  in  either  an  elastic  or  a 
dissipative  region,  removed  from  an  interface  or  discontin¬ 
uity.  Figure  11  illustrates  a  step  from  t  to  t  +  At  . 

The  values  of  all  functions  are  known  for  t  =  t  ,  i.e., 

1 
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along  A-E.  At  point  P  ^  there  are  five  unknowns,  which  may 
be  obtained  from  the  respective  five  equations  along  the 
characteristics,  subject  to  check  on  the  yield  condition,  or 
on  A  )>  0  ,  as  appropriate.  There  is  no  difficulty  at  a 
point  where  yield  has  not  been  reached,  the  elastic  relations 
must  certainly  apply.  However,  disregarding  discontinuities, 
consider  a  situation  where  the  yield  condition  is  satisfied 
for  a  portion  of  the  range  A-E,  indicating  that  there  might 
be  an  elastic-plastic  boundary.  To  locate  a  point  P  on  the 
boundary  and  the  inclination  a  of  the  latter.  Fig.  12-a, 
one  can  combine  the  conditions  along  those  elastic  and  plas¬ 
tic  characteristics  which  lie  in  the  proper  regions,  but, 
depending  on  the  value  of  Qj  there  are  a  variety  of  possi¬ 
bilities.  There  are  six  unknown  quantities,  the  inclination 
a  ,  three  stresses,  and  two  velocities,  while  the  yield  con¬ 
dition  and  five  or  six  characteristic  conditions  are  avail¬ 
able.  Over-determination,  Fig.  12-b,  in  case  the  seven  con¬ 
ditions  apply  is  no  difficulty.  It  may  simply  mean  that  the 
particular  location  may  not  occur.  There  are,  however, 
values  of  a  as  shown  in  Pig.  12-c  where  six  conditions 
apply,  one  of  which  is  an  inequality,  with  a  possibility  of 
indeterminancy.  Before  proceeding  with  a  numerical  analysis, 
it  is  necessary  to  ascertain  the  existence  and  uniqueness 
of  the  solutions  for  the  location  of,  and  stresses  at  a 
boundary  point  P  from  the  conditions  available  for  all  pos¬ 
sible  situations.  If  this  can  be  achieved,  the  additional 
possibility  of  discontinuities,  the  velocities  of  propaga¬ 
tion  of  which  are  known,  should  be  no  serious  obstacle. 
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AND  -O',  ^  -Cro<  -cr,  FOR  07  ,  O',  ,  SEE  EQ  S. ( 4 -4)  ,  ( 4 -  7  ) 


54 


NOTE  :  MATERIAL 


TABLE 
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NOTE-,  material  BECOMES  PLASTIC  AT  P-FRONT 
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FIG.  I  half-space  subjected  TO  UNIFORM  LOADS 
O*o(t)  AND  To(t)  AT  THE  SURFACE  X=0 


-57- 


i.5 


(no  upper  limit  on 


REGION  FOR  CASE  A- 3 
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YIELD  AT  P-FRONT 
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NO  SOLUTION 
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DEPENDENCE  OF  TYPE  OF  SOLUTION  ON  SURFACE  STRESSES 
OTo  ,  To  WHEN  (3  >5  (DRAWN  FOR  (3  =  5) 
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TYPICAL  STRESS  HISTORY  FOR  CASE  A. 2 


FIG.  7  DEPENDENDE  OF  TYPE  OF  SOLUTION  ON  SURFACE  STRESSES 
03  ,  to  WHEN  /3  <  3  (DRAWN  FOR  /3  =  2  ) 
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FIG.  9  half-space  loaded  BY  MOVING 
STEP  PRESSURE  ON  SURFACE 
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Department  of  the  Navy 
War-hington  25,  D.  C, 

Attn :  Code  32? 

Code  348 
Code  371 
Code  420 
Code  423 
Code  442 
Code  421 

Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D,  C. 

Attn:  AD-2 
ad-22 
RS-7 
TD-42 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 


Attn:  Ad3  (l) 

Re  ( 1 1 

Re3  (1) 

Re5  (1) 

ReN  ( 1 ) 


ili 


Commanding  Officer 
Watertown  Arsenal 
Watertown,  Massachusetts 
Attn;  Laboratory  Div.  (l) 


Director  of  Naval  Intelligence 
Navy  Department 
Washington  25,  D.  C. 

Attn:  Op-922V  (1) 


Commanding  Officer 

Frankford  Arsenal 

Cridesburg  Station 

Philadelphia  37,  Pa- 

Attn:  Leiboratory  Div.  (1) 

Office  o*"  Ordnance  Research 

2127  Myrtle  Drive 

Duke  Station 

Durham,  North  Carolira 

Attn:  Div.  of  Engin.  Science8(l) 


Chief,  Bureau  of  Yards  eind  Docks 
Department  of  the  Navy 
Washington  25,  D,  C. 

Attn;  Code  D-213  (1) 

Code  D-222  (l) 

Code  D-41CX:  (1) 

Code  D-440 

Code  D-500  (1) 


Commanding  Officer 

Engin.  Research  Develop.  Lab. 

Fort  Belvoir,  Virginia  (1) 


NAVY  (Continued) 


Commanding  Officer  and  Director 
David  Taylor  Model  Basin 
Washington  7.  D.  C. 


Attn:  Code  140  (1 

Code  600  ( 1 

Code  700  (1 

Code  720  ( 1 

Code  726  ( 1 

Code  731  (l 

Code  740  (2 


Commander 

U.  S.  Naval  Ordnance  Laboratory 


White  Oak,  Maryland 
Attn:  Tech.  Library 

Tech.  Evaluation  Dep.  (1) 
EE  (l) 
EH  (  1 } 
R  (1) 


Director 

Material  Laboratory 
New  York  Naval  Shipyard 
Brooklyn  1,  New  York  (1) 


Commanding  Officer  and  Director 
U.  S.  Nava!  Electronics  Lab. 

San  Diego  62,  California 
Attn:  Code  4223 


Of f icer-in-Charge 
Naval  Civil  Engin.  Research 
Naval  Civil  Engin.  Research  .and 
Evaluation  Laboratory 
U.  S.  Naval  Construction 
Battalion  Center 
Port  Hueneme,  California 

Attn:  Code  763 


Director 

Naval  Air  Experimental  Station 
Naval  Air  Material  Center 
Naval  Base 
Philadelphia  12,  Pa 
Attn:  Materials  Laboratory  (1) 
Structures  LaJa.  (l) 


Of freer- in-Charge 

Underwater  Explosion  1. -search  Div. 
Norfolk  Naval  Shipyard 
Portsmouth,  Virginal 
Attn:  Dr.  A.H.  Keil  (2) 


Commander 

U.S.  Naval  Providing  Grounds 
Dahlgren,  Virginia 

Superintendent 
Naval  Gun  Factory 
Washington  26,  D.  C. 

Commander 

Naval  Ordnance  Test  Station 
Inyokern,  China  Lake,  California 
Attn:  Physics  Division 

Mechanics  Branch 

Commander 

Naval  Ordnance  Test  Station 
Underwater  Ordnance  Division 
32C2  E.  Foothill  Boulevard 
Pasadena  6,  California 
Attn:  Structures  Division 

Commanding  Officer  and  Director 
Naval  Engineering  Experiment  Station 
Annapolis,  Maryland  (1) 

Super intendent 

Naval  Post  Graduate  School 

Monterey,  California  (1) 

Commandant 

Marine  Corps  Schools 

Quantico,  Virginia 

Attn:  Director,  Marine  Corps 

Development  Center  (1) 

AIR  FORCE 

Commanding  General 
U.  .S.  .^ir  Force 
Washington  25,  D.  C. 

Attn:  Research  and  Development  Div.{l) 

Commander 

Air  Materiel  Command 

Wright- Pat ter son  Air  Force  Base 

Dayton,  Ohio 

Attn:  WCOSI  (l) 

Commanac  r 

U.  S.  A.”  force  Institute  of  Technology 
Wright- Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  Chief,  Applied  Mech.  Group  (l) 

Director  of  Intelligence 
Headquarters,  U.  S.  Air  Force 
Washington  25,  D.  C. 

Attn:  P,  V.  Br.  (Air  Targets  Div.)  fl) 

AF0IN-1B2  (2) 


(1'/ 

(1) 


Commander 

Air  Research  and  Development 
Command 

P.  O,  Box  1395 

Attn:  RDMPE  (l) 

Commander 

W.VDD 

Wright-Patterjon  Air  Force  Base 
Dayton j  Ohio 

Attn:  WWRC  (1) 

WWRMDS  ( • \ 

WWRMDD  (1) 

Commanding  Officer 
USNNOEU 

Kirtland  Air  Force  Base 
Albuquerque,  New  Mexico 
Attn ;  Code  20 

(Dr.  J.  N.  Brennan)  (l) 

OTHER  GOVERNMENT  ACTIVITIES 

U.  S.  Atomic  Energy  Commission 
VJashington  25,  D.  C. 

Attn:  Director  of  Research  (2) 

Director 

National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Division  of  Mechanics  (l) 
Engin.  Mechanics  Sec.  (11 
Aircraft  Structures  (l) 

Commander 

U.  S.  Coast  Guard 

1300  E.  Street,  N.  W. 

Washington  25,  D.  C. 

Attn:  Chief,  Testing  and 

Developmenc  Div.  (l) 

U.S.  Maritime  Administration 
General  Administration  Office  Bldg. 
Washington  25,  D.  C. 

Attn:  Chief,  Div.  of  Preliminary 
Design  (ll 

National  Advisory  Committee  for 
Aeronautics 
1512  H.  Street,  N.W. 

Wdshing  Loll  25,  D.  C. 

Attn:  Loads  and  Structures  Div. 

(2) 

Director 

Langley  Aeronautical  Laboratory 

Langley  Field,  Virginia 

Attn:  Structures  Div.  (2) 


Director 

Forest  Products  Laboratory 

Madison,  Wisconsin  (l) 

Civil  Aeronautical  Administration 
Denar tmenL  of  Commerce 
Washington  25,  D.  C, 

Attn:  Chief,  Airframes  and  Equip. Br.(l) 

National  Sciences  Foundation 
1520  H.  Street.  N.  W. 

Washington  25,  D.  C. 

Attn:  Engin.  Science  Division  (l) 

National  Academy  of  Science 
2101  Constitution  Avenue 
Washington  25,  D.  C. 

Attn:  Tech.  Director,  Committee  On 

Ships’  Structural  Design  (l) 

Executive  Secretary,  Committee 
on  Underwater  Warfare  (l) 

Director,  Operations  Research  Office 

John  Hopkins  University 

7100  Connecticut  Avenue 

Chevy  Chase,  Maryland 

Washington  15,  D.  c.  (1) 

Dr.  Alvin  c.  Graves,  Director 
J-Division,  Los  Alamos  Scientific  Lab. 

P.  O,  Box  166 j 

Los  Alamos,  New  Mexico  (1) 

U.  S.  Atomic  Energy  Commision 
Classified  Technical  Library 
1901  Constitution  Avenue,  N.  W. 
Washington,  D.  C. 

Attn:  Mrs.  Jean  M.  O'Leary 

for  Dr.  Paul  C.  Pine  (l) 

U.S.  Atomic  Energy  Commission 
Classified  Tech.  Library 
Tech.  Information  Service 
1901  Constitution  Avenue,  N.W. 

tn:  Mrs.  Jean  M. O'Leary  (l) 

Sandia  Corporation,  Sandia  Base 

Albuquerque,  New  Mexico 

Attn:  Dr.  Walter  A.  MacNair  (l) 

LegisidLive  Reference  Service 
Library  of  Congress 
Washington  25,  D.  C. 

Attn:  Dr.  E.  Wenk  (l) 

U.S.  Atomic  Energy  Commission 
1901  Constitution  Avenue,  N.W. 
Washington,  D.  C. 


OTHER  GOVERNMENT  ACTIVITIES  (Continued) 


Massachusetts  Institute  of 
Technology 

Cainbridge  39>  Massachusetts 
Attn:  Dr,  Charles  H.  Norris  (l) 

INVESTIGATORS  ACTIVELY  ENGAGED 
IN  RELATED  RESEARCH _ 


Professor  Lynn  S.  Beedle 
Fritz  Engineering  Laboratory 
Lehigh  University 

Bethlehem,  Pennsylvania  (l) 

Professor  H.  H.  Bleich 
Department  of  Civil  Engineering 
Columbia  University 
New  York  2?)  New  York  (1) 

Professor  R.  L.  BispJinghoff 
Dept,  of  Aeronautical  Engineering 
Massachusetts  Institute  of  Techno. 

Cambridge  39>  Massachusetts  (l) 

Professor  B.  A.  Boley 
Department  of  Civil  Engineering 
Columbia  University 
New  York  27,  N.Y. 

Professor  Eugene  J.  Brunei le.  Jr. 
Department  of  Aeronautical  Engin. 
Princeton  University 
Princeton,  New  Jersey  (1 

Professor  G.F.  Carrier 

Pierce  Hall 

Harvard  University 

Cambridge  3®>  Massachusetts  (  1 

Professor  J.E.  Cermak 
Department  of  Civil  Engin. 

Colorado  state  University 

Fort  Collins,  Colorado  (1 

Professor  Herbert  Deresiewicz 
Department  of  C  vil  Engin. 

Columbia  University 

632  West  125th  street 

New  York  27,  N.Y.  (!' 

Professor  Lloyd  Donnell 
Department  of  Mechanics 
Illinois  Institute  of  Technology 

Chicago  I6,  Illinois  ( 1 


Professor  D.  C,  Drucker,  Chairman 
Division  of  Engineering 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  A.  C.  Eringen 

Department  of  Aeronautical  Engineering 

Purdue  University 

Lafayette,  Indiana  (i) 

Professor  W.  Flugge 

Department  of  Mechanical  Engineering 
Stanford  University 

Stanford,  California  (1) 

Mr.  Martin  Goland 

Midwest  Research  Institute 

40^9  Pennsylvania  Avenue 

Kansas  City  2,  Missouri  (l) 

Professor  J.  N.  Goodier 

Department  of  Mechanical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Professor  L.E.  Goodman 
Engineer  Experiment  station 
University  of  Minnesota 
Minneapolis,  Minnesota  (1^ 

Professor  W.J.  Hall 

Department  of  Civil  Engineering 

University  of  Illinois 

Urbana,  Illinois  (l) 

Professor  R.P.  Harrington,  Head 
Department  of  Aeronautical 
Engineering 

University  of  Cincinnati 

Cincinnati  21,  Ohio  (l) 

Professor  M.  Hetenyi 

The  Technological  Instil  ate 

Northwestern  University 

Evanston,  Illinois  (1) 

Professor  P.G.  Hodge 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Technology  Center 

Chicago  16,  Illinois  (l) 

^JfOfGSSOJT  N  .  J  •  Ho ^ f 
Stanford  University 
Stanford,  California 


(1) 


INVESTIGATORS  ACTIVELY  ENGAGED 
IN  RELATED  RESEARCH  (Continued) 


Professor  W.  H.  Koppmann,  II 
Department  of  Mechanical  Engr. 
John  HopKins  University 
Baltimore,  Maryland  (1) 

Professor  Bruce  G.  Johnson 

University  of  Michigan 

Ann  Arbor,  Michigan  (1) 

Professor  J,  Kempner 
Department  of  Aeronautical  Engr. 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
99  Livingston  Street 
Brooklyn  2,  New  York  (1) 

Professor  H.  L.  Langhaar 
Department  of  Theoretical 
and  Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  B.  J.  Lazan,  Director 
Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis  14,  Minnesota  (1) 

Professor  E.  H.  Lee 
Stanford  University 
Stanford.  California 

(1) 

Professor  George  H.  Lee 
Director  of  Research 
Rensselaer  Polytechnic  Institute 
Troy ,  New  York 

Mr.  M.  M.  Lemcoe 

Southwest  Research  Institute 

8500  Culebra  Road 

San  Antonio  6,  Texas  (1) 

Professor  Paul  Lieber 
Geology  Department 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (l) 

Professor  Hsu  Lo 

School  of  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Professor  R.  D.  Mindlin 
Department  of  Civil  Engineering 
Columbia  University 
New  York  27,  N.  Y.  (1) 


Dr.  A.  Nadai 

136  Cherry  Valley  Road 

Pittsburgh  21,  Pennsylvania  (l) 

Professor  William  A.  Nash 
Department  of  Engineering  Mechanics 
University  of  Florida 

Gainesville,  Florida  (l) 

Professor  N.  M.  Newmark 

Department  of  Civil  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Professor  Aris  Phillisps 
Department  of  Civil  Engineering 
15  Prospect  Street 
Yale  University 

New  Haven,  Connecticut  (l) 

Professor  W.  Prager,  Chairman 
Physical  Sciences  Council 
Brown  University 

Providence  12,  Rhode  Island  (l) 

Professor  E.  Reissner 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (1) 

Professor  M.  A.  Sadowsky 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 

Dr,  B.  w.  Shaffer 

Department  of  Mechanical  Engineering 
New  York  University 
45  Fourth  Avenue 

New  York  5>3>  New  York  (l) 

Professor  C.  B.  Smith 
College  of  Arts  and  Science 
Department  of  Mathematics 
Walker  Hall 
University  of  Florida 

Gainesville,  Florida  (l) 

Professor  J,  Stallmeyer 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Eli  Sternberg 
Stanford  University 
Stanford  California 


(1) 


IN'/ESTIGATORS  ACTI^/ELY  ENGAGED 
IN  RELATED  RESEARCH  ^Contiiiued  ’ 

Professor  S.P.  Timoshenko 
School  of  Engineering 
Stanford  University 
Stanford,  California  , 1  ; 

Professor  A. A.  Velestos 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  vl 

Professor  Enrico  Volterra 
Department  of  Engineering  Mech . 
University  of  Texas 
Austin  12,  Texas  ,1' 

Dr.  Dana  Young 

Southwest  Research  Institute 
San  Antonio.  Texas  1 

Project  Staff  liv 


